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Abstract 

Wightman function, the vacuum expectation values of the field square and the energy- 
momentum tensor are evaluated for a scalar field obeying the Robin boundary conditions on 
two spherical branes in (D + l)-dimcnsional Rindler-like spacetime Ri x S ^ 1 , with a two- 
dimensional Rindlcr spacetime Ri. This spacetime approximates the near horizon geometry 
of (D + l)-dimcnsional black hole. By using the generalized Abel-Plana formula, the vacuum 
expectation values are presented as the sum of single brane and second brane induced parts. 
Various limiting cases are studied. The vacuum forces acting on the branes are decomposed 
into the self-action and interaction terms. The interaction forces are investigated as functions 
of the brane locations and coefficients in the boundary conditions. 

PACS number(s): 03.70.+k, 04.62,+v, ll.10.Kk 

1 Introduction 

Motivated by string/M theory, the AdS/CFT correspondence, and the hierarchy problem of 
particle physics, braneworld models were studied actively in recent years (for a review see pQ). 
In these models, our world is represented by a sub-manifold, a three-brane, embedded in a 
higher dimensional spacetime. In particular, a well studied example is when the bulk is an 
AdS space. The braneworld corresponds to a manifold with boundaries and all fields which 
propagate in the bulk will give Casimir-type contributions to the vacuum energy, and as a result 
to the vacuum forces acting on the branes. In dependence of the type of a field and boundary 
conditions imposed, these forces can either stabilize or destabilize the braneworld. In addition, 
the Casimir energy gives a contribution to both the brane and bulk cosmological constants 
and, hence, has to be taken into account in the self-consistent formulation of the braneworld 
dynamics. Motivated by these, the role of quantum effects in braneworld scenarios has received 
a great deal of attention. For a conformally coupled scalar this effect was initially studied in 
Ref. [2] in the context of M-theory, and subsequently in Refs. |Hj for a background Randall- 
Sundrum geometry. The models with dS and AdS branes, and higher dimensional brane models 
are considered as well [I]. 

* E- mail : sahary an@ server . phy sdep . r . am 
^E-mail: rezakord@ipm.ir 
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In view of these recent developments, it seems interesting to generalize the study of quantum 
effects to other types of bulk spacetimes. In particular, it is of interest to consider non-Poincare 
invariant braneworlds, both to better understand the mechanism of localized gravity and for 
possible cosmological applications. Bulk geometries generated by higher-dimensional black holes 
are of special interest. In these models, the tension and the position of the brane are tuned in 
terms of black hole mass and cosmological constant and brane gravity trapping occurs in just the 
same way as in the Randall-Sundrum model. Braneworlds in the background of the AdS black 
hole were studied in J5]. Like pure AdS space the AdS black hole may be superstring vacuum. 
It is of interest to note that the phase transitions which may be interpreted as confinement- 
deconfinement transition in AdS/CFT setup may occur between pure AdS and AdS black hole 
0. Though, in the generic black hole background the investigation of brane-induced quantum 
effects is technically complicated, the exact analytical results can be obtained in the near horizon 
and large mass limit when the brane is close to the black hole horizon. In this limit the black 
hole geometry may be approximated by the Rindler-like manifold (for some investigations of 
quantum effects on background of Rindler-like spacetimes see [7] and references therein). 

In the previous papers [El El we have considered the vacuum densities induced by a spherical 
brane in the bulk Ri x S D_1 , where Ri is a two-dimensional Rindler spacetime. Continuing 
in this direction, in the present paper we investigate the Wightman function, the vacuum ex- 
pectation values of the field square and the energy-momentum tensor for a scalar field with an 
arbitrary curvature coupling parameter for two spherical branes on the same bulk. Though the 
corresponding operators are local, due to the global nature of the vacuum, these expectation 
values describe the global properties of the bulk and carry an important information about the 
physical structure of the quantum field at a given point. The expectation value of the energy- 
momentum tensor acts as the source of gravity in the Einstein equations and, hence, plays 
an important role in modelling a self-consistent dynamics involving the gravitational field. In 
addition to applications in braneworld models on the AdS black hole bulk, the problem under 
consideration is also of separate interest as an example with gravitational and boundary-induced 
polarizations of the vacuum, where all calculations can be performed in a closed form. Note 
that the vacuum densities induced by a single and two parallel flat branes in the bulk geometry 
Ri x R D_1 for both scalar and electromagnetic fields are investigated in jlUl lllj . 

The paper is organized as follows. In the next section we consider the positive frequency 
Wightman functions in the region between two branes. On the basis of the generalized Abel- 
Plana formula, we present this function in the form of the sum of single brane and second brane 
induced parts. By using expression for the Wightman function, in section El we investigate 
the vacuum expectation values of the field square and the energy-momentum tensor. Various 
limiting cases of the general formulae are studied. In section 0] the vacuum forces acting on the 
branes due to the presence of the second brane are evaluated by making use of the expression 
for the radial vacuum stress. The main results of the paper are summarized in section 

2 Wightman function 

We consider a scalar field ip(x) propagating on background of (D + l)-dimensional Rindler-like 
spacetime Ri x S D_1 . The corresponding metric is described by the line element 

ds 2 = fdr 2 -df -r\dY? D _ x , (1) 

with the Rindler-like (r, £) part and dY? D _ x is the line element for the space with positive constant 
curvature with the Ricci scalar R = (D — 2)(D — \)jr\. Line element Q describes the near 
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horizon geometry of (D + l)-dimensional topological black hole with the line element 12 



ds 2 = A H {r)dt 2 - -^L— - r 2 dS 2 D ^, (2) 

where Ah(t) = k + r 2 /l 2 — /l 2 r n , n = D — 2, and the parameter k classifies the horizon 
topology, with k = 0, —1, 1 corresponding to flat, hyperbolic, and elliptic horizons, respectively. 
The parameter I is related to the bulk cosmological constant and the parameter r$ depends on 
the mass of the black hole. In the non extremal case the function Au(r) has a simple zero at 
r = rjj, and in the near horizon limit, introducing new coordinates r and p in accordance with 

t = A' H (r H )t/2, r-r H = A' H {r H )e/4, (3) 

the line element is written in the form l|T|). Note that for a (D + l)-dimensional Schwarzschild 
black hole one has An(r) = 1 — (rn/r) n and, hence, A' h (th) = n/rn- 
The field equation is in the form 

(g ik ViV fe + m 2 + CR) <p(x) = 0, (4) 

where £ is the curvature coupling parameter. Below we will assume that the field satisfies the 
Robin boundary conditions on the hypersurfaces £ = a and £ = b, a < b, 



0_ 



= 0, j = a, b, (5) 



with constant coefficients Aj and Bj. The Dirichlet and Neumann boundary conditions are 
obtained as special cases. In accordance with (j3J), the hypersurface £ = j corresponds to the 
spherical shell near the black hole horizon with the radius rj =ru + A' H (rH)j 2 /4. 

The branes divide the bulk into three regions corresponding to < £ < a, a < £ < 6, and 
b < £ < oo. In general, the coefficients in the boundary conditions © can be different for 
separate regions. In the corresponding braneworld scenario based on the orbifolded version of 
the model the region between the branes is employed only and the ratio Aj/ Bj for untwisted 
bulk scalars is related to the brane mass parameters Cj of the field by the formula [H] 

A i - 1 f, c 



Bj 2{ Cj -y^ j=a ' k (6) 

For a twisted scalar the Dirichlet boundary conditions are obtained on both branes. 

To evaluate the vacuum expectation values (VEVs) of the field square and the energy- 
momentum tensor we need a complete set of eigenfunctions satisfying the boundary conditions 
In accordance with the problem symmetry, below we shall use the hyper spherical angular 
coordinates 4>) = (61,62, ... , 6 n , (j)) on S ' 1 with ^ 0^ ^ 7r, k = 1, . . . , n, and 0^0^ 27T. 
In these coordinates the eigenfunctions in the region between the branes can be written in the 
form 

<^(x) = C a zS(A,e,A,6)r(m fc ;i?,0)e- <WT , (7) 

where rrik = (mo = I, mi, . . . m n ), and mi,m2, ■ ■ ■ m n are integers such that ^ m n -\ ^ • • • ^ 
m\ ^ I, — m n _i ^ m n ^ m n _i. In Eq. Q Y(mk;$, <fi) is the spherical harmonic of degree I [Tlj . 
and 

Z$(u,v) = lS(v)KUu) - K$(v)I i0J (u), j = a,b, (8) 
with Iiuj(x) and Ki UJ (x) being the modified Bessel functions with the imaginary order, 

Xl = —Jl(l+n)+Cn(n+l)+m 2 r 2 H . (9) 
r H v 
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Here and below for a given function f(z) we use the barred notations 



fO)(z)=A J f(z) + %f'(z), j = a,b. (10) 



Functions (J7J) satisfy the boundary condition on the brane £ = b. From the boundary 
condition on the brane £ = a we find that the possible values for u are roots to the equation 

Z iw (X l a,X l b) = 0, (11) 

with the notation 

Z w (u,v)= (v)K^ («) - K-W («) («) . (12) 

For a fixed A;, the equation Hll|) has an infinite set of real solutions with respect to uj. We will 
denote them by to n = u n (Xia, Xfi), u n > 0, n = 1, 2, . . ., and will assume that they are arranged 
in the ascending order uj n < u> n+ \. In addition to the real zeros, in dependence of the values 
of the ratios jAj/Bj, equation (Jll[) can have a finite set of purely imaginary solutions. The 
presence of such solutions leads to the modes with an imaginary frequency and, hence, to the 
unstable vacuum. In the consideration below we will assume the values of the coefficients in 
boundary conditions © for which the imaginary solutions are absent and the vacuum is stable. 
The coefficient C a in can be found from the normalization condition 

/rb j y 

dTt J —tp a d T cp* a , = i5 aa t. (13) 

where the integration goes over the region between two spheres. Substituting eigenfunctions l(7|) 
and using the relation j \Y(mk] (j))\ 2 dVl = N{m,k) for spherical harmonics, one finds 



C 2 _ r]f D I^(X ia ) 



N(m k )lZ'(Xib)£Z iuJ ^ia,Xib) 



(14) 



The explicit form for N{rrik) is given in Jl] and will not be necessary for the following consid- 
erations in this paper. 

First of all we evaluate the positive frequency Wightman function 

G+(x,x') = (0\(p{x)(p(x')\0), (15) 

where |0) is the amplitude for the corresponding vacuum state. The VEVs for the field square 
and the energy-momentum tensor are obtained from this function in the coincidence limit of the 
arguments. In addition, the Wightman function determines the response of a particle detector 
in given state of motion. By expanding the field operator over eigenfunctions and using the 
commutation relations one can see that 

G + {x,x') = Y,^{x)^ a {x'). (16) 

a 

Substituting eigenfunctions into this mode sum formula and by making use of the addition 
theorem 

N(m k ) nmk >*^ ) = ^ C > {COSO) > (17) 
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for the Wightman function in the region between the branes one finds 

x^A^A^^A^A^I . (18) 

\(jJ=U} n 

In this formula, Sd = 2ir D l 2 /F(D/2) is the volume of the unit (D — l)-sphere, C? (x) is the 
Gegenbauer polynomial of degree I and order n/2, 9 is the angle between directions (1?, (f>) and 

As the normal modes uj n are not explicitly known and the terms with large n are highly 
oscillatory, the Wightman function in the form is not convenient. For the further evaluation 
we apply to the sum over n the summation formula derived in Ref. on the basis of the 
generalized Abel-Plana formula |15| : 

ti&ZU","^, Jo » 2 Jo 2*Z z (u,v) •*■>-*■>'> I 

As a function F(z) in this formula we choose 

The conditions for the formula (|19|) to be valid are satisfied if a 2 e' T_r ' < For the Wightman 
function one obtains the expression 

= G + (x,x';6) - - JL - r - V (2/ + n)C, n/2 (cos 0) / QUXia, \ t b) 

xZ< 6 )(A,£, A,6)ZW(A,^, A/6) cosher - r')], (21) 

7 (a) (ii) 

£M^) = (fe) — - — -• (22) 

„1-D 00 „oo 

G+(a;,x / ;&) = J 1 - V(2Z + n)C, n/2 (cos 9) / dwsinh(7ra;) 
7r 2 n5 D ^ Jo 



where 



In Eq. (J2U) 



^(^6)111(^6) 

is the Wightman function in the region £ < 6 for a single brane at £ = 6 and the second term 
on the right is induced by the presence of the brane at £ = a. The function (|2*Tj) is investigated 
in Ref. [8] and can be presented in the form 

G+(x,x';b) = G+(x,x') + (<p(x)^(x'))( b \ (24) 

where Gq(x,x') is the Wightman function for the geometry without boundaries and the part 

M*M*')> (6) = -^jr(2l + n)cf 2 (cos9) r^^M 
™ S v\To J ° IiW) 

cosher - r')] (25) 
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is induced in the region £ < b by the presence of the brane at £ = b. Note that the representation 
(|24|) with (|25jl is valid under the assumption ££' < 6 2 e' T_T I As it has been shown in jS], the 
Wightman function for the boundary-free geometry may be written in the form 

G+(x,x') = sO-HFc- E( 2Z + n K /2 ( cos 

^ n D ;=o 

/>oo 

x / dwe"^ cos[o;(T-r / )]^(A / e)^(A i O, (26) 
J o 

where G[j~(x, x') is the Wightman function for the bulk geometry R 2 x Outside the horizon 

the divergences in the coincidence limit of the expression on the right of (|26j) are contained in 
the first term. 

It can be seen that the Wightman function in the region between the branes can be also 
presented in the form 

r l-D oo roc 

G + (x,x') = G + (x,x';a) - -Z— V (2/ + n)C, n/2 (cosfl) / duj n a J\ ia , Xib) 

x^(¥, \ ia )Z^(X^, A/a) cosher - r% (27) 



with the notation 



In this representation, 



^(u,v)= u U . (28) 

{U, V ) 



G + (x, x'; a) = G+(x, x') + {<p(x)<p(x'))M (29) 
is the Wightman function in the region £ > a for a single brane at £ = a, and 



r l-D OO ,00 f( a V\ n\ 

(^» (,) = -^vE( 2/ +^r /2 (cos^) / 

^ «"i o) (A,a) 



x^(A/0^(A^') cosh[u,(r - r')]. (30) 

Two representations of the Wightman function, given by Eqs. (|21j) and (|27|) . are obtained from 
each other by the replacements 

a +± 6, J w (31) 

In the coincidence limit the second term on the right of formula (|21j) is finite on the brane £ = b 
and diverges on the brane at £ = a, whereas the second term on the right of Eq. (|27|) is finite on 
the brane £ = a and is divergent for £ = b. Consequently, the forms (|21|) and (|27|) are convenient 
for the investigations of the VEVs near the branes £, = b and (, = a, respectively. 

We have investigated the Whightman function in the region between two branes for an 
arbitrary ratio of boundary coefficients Aj/Bj. Note that in the orbifolded version of the model 
the integration in the normalization integral goes over two copies of the bulk manifold. This 
leads to the additional coefficient 1/2 in the expression 1)14(1 for the normalization coefficient C a . 
Hence, the Whightman function in the orbifolded braneworld case is given by formula 1)2 ljl with 
an additional factor 1/2 in the second term on the right and in formula (|23fl . As it has been 
mentioned above this function corresponds to the braneworld in the AdS black hole bulk in the 
limit when the branes are close to the black hole horizon. 
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3 Casimir densities 



3.1 VEV for the field square 

In this section we will consider the VEVs of the field square and the energy-momentum tensor 
in the region between the branes. In the coincidence limit, taking into account the relation 
C?V 2 (1) _ _|_ n yp(n)l\, from the formulae for the Wightman function one obtains two 
equivalent forms for the VEV of the field square: 

(0\<f?\0) = (0 |^|Oo) + (^) (j) 

— / dun^X^X^Z^iX^Xd), (32) 

corresponding to j = a and j = 6, and |Oo) is the amplitude for the vacuum without boundaries, 

D '=^ +D -v T -m^w (33) 

is the degeneracy of each angular mode with given /. The VEV (0 |^ 2 |0 ) is obtained from 
the corresponding Wightman function given by (|26|), For the points outside the horizon, the 
renormalization procedure is needed for the first term on the right only, which corresponds to 
the VEV in the geometry R 2 x S " 1 . This procedure is realized in |Sj on the base of the zeta 
function technique. 

In Eq. (|32[). the part {(f 2 )^ is induced by a single brane at £ = j when the second brane is 
absent. For the geometry of a single brane at £ = a, from Q3UJI one has 

W'^-^EA (34) 
ttSd Jo K^ ] {X ia ) 



1=0 



The expression for {(f 2 )^ is obtained from (|34jl by replacements 1)31(1 and is investigated in |Sj. 
The last term on the right of formula 1)32)1 is induced by the presence of the second brane. It 
is finite on the brane at £ = j and diverges for the points on the other brane. By taking into 
account the relation zjj\u,u) = Bj/j, we see that for the Dirichlet boundary condition this 
term vanishes on the brane £ = j. 

Let us consider the behavior of the single brane part Q34j) in asymptotic regions of the 
parameters. In the limit £ — > a this part diverges and, hence, for points near the brane the main 
contribution comes from large values I. By making use of the corresponding uniform asymptotic 
expansions for the modified Bessel functions, to the leading order we find 

y>W « - X ( ^ n - , (35) 
(4vr)— (t-a) - 1 

where <5e a = 1 for B a = and <5e a = —1 for B a ^ 0. Hence, near the brane the brane- induced 
part is negative for the Dirichlet boundary condition and is positive for non-Dirichlet boundary 
condition. At large distances from the brane, £ 3> rjj, the dominant contribution into (|34j) 
comes from the I = term and in the leading order we have 

(^-if^r^l^. (36 ) 

2Sz?A £ Jo K^ a) (Xoa) 

In the limit when the position of the brane tends to the horizon, a — > 0, with fixed £, we use 
the formulae for the modified Bessel functions with small values of the argument. The main 
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contribution into the integral comes from the lower limit of the integration and we obtain the 
formula 



2\(o) 



£a^o 2 (M)- 



(37) 



27r5 D ln (r H /a) £^ 

In the limit r — ► the curvature for the background spacetime is large. In this limit A; 
is also large. The exception is the term I = for a minimally coupled scalar field for which 
Ao = m. For large values A; the main contribution into the integral in (|34|) comes from large 
values u. Introducing a new integration variable x = w/A/a, we estimate the integral by the 
Laplace method. This leads to the following result 



2\(o) 



SB a r)T D a e W [-2VCn(n + m ~ a)/r H ] 



(38) 



For a minimally coupled scalar field the contribution of the terms with I ^ 1 is suppressed by 
the factor e -2A '^~ a ) and the dominant contribution comes from the I = term: 



2\(tt) 



' H 

ttSd 



duj- 



iffi (ma) 



K 



(39) 



ma l 



As we see, the behavior of the VEV in the high curvature regime is essentially different for 
minimally and non-minimally coupled fields. 

In the near-horizon limit, a, £ <C ru, the main contribution into the sum over I in (|34j) comes 
from large values Z corresponding to Z < th/H — o). To the leading order we can replace the 
summation over I by the integration in accordance with the formula 



oo 

E 

1=0 



Dif(Xi) 



2r 



D-l 
H 



T(D - 1) 



dk k D - 2 f(Vk 2 + m 2 ). 



(40) 



Now it is easily seen that from (|34j) we obtain the corresponding result for the plate uniformly 
accelerated through the Fulling-Rindler vacuum. 

In the geometry of two branes, extracting the contribution from the second brane, we can 
write the expression (|32|) for the VEV in the symmetric form 



(0|^ 2 |0) = (Oo|v 2 |Oo)+ 2J(f? 

j=a,b 



j2 (ip 2 )^ + (ip 2 ^" 1 " 



(41) 



with the interference part 



2\(ab) 



1-D oo 
H 



ttS d 



i=o Jo 



lL b \\ib)Z w (Xia,Xib) 



K^\x ia ) 



(42) 



An equivalent form for this part is obtained with the replacements (|31|) in the integrand. The 
interference term ()42l) is finite for all values of £ in the range a ^ £ ^ b, including the points on 
the branes. The surface divergences are contained in the single brane parts only. 

Let us consider the behavior of the interference part in the VEV of the field square in limiting 
regions for values of the parameters. First of all, it can be seen that in the limit a — > b, to the 
leading order the result for the parallel plates in the Minkowski bulk is obtained. When the left 
brane tends to the horizon, a — > 0, the dominant contribution comes from the lower limit of the 
integration in (j42fl . and we have 



2\ (ab) 



27rS , £)ln (rn/a) 



(b) 



(Aj&) 



1=0 



2K (A,0 



K, 



(6) 



(Mb) 
u 



I { o h \kb) 



k{XiO- (43) 
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In the limit b — ► oo for fixed values a and £, the main contribution comes from the lowest mode 
1 = 0, and to the leading order one finds 



2\(a6) 



S^g" 1 Aj, + B b X Jo 



00 , /i a) (A a) 
aw- 



^ a) (A a) 



Zi a) (A a) 
i^ a) (A a) 



(44) 

with the exponentially suppressed interference part. The behavior of the interference part in 
the limit r# — * can be investigated in the way similar to that for a single brane part. For a 
non-minimally coupled scalar field the interference part is dominated by the I = term and is 
suppressed by the factor exp[— 2^/Qn{n + l)(b — a)/V#]. For a minimally coupled scalar field 
the leading term is given by the I = summand with Ao = m and the interference part behaves 
as r]f D . In the near-horizon limit, a,b <C rn, replacing the summation by the integration in 
accordance with formula Q4UJI . it can be seen that from (|42[l the result for the geometry of two 
parallel plates uniformly accelerated through the Fulling-Rindler vacuum is obtained. 



3.2 Energy-momentum tensor 

The VEV of the energy-momentum tensor is expressed in terms of the Wightman function as 



(0|T«|0) = lim did' k G + (x,x') + 



(0|^|0), (45) 



where Rik is the Ricci tensor for the bulk geometry. Making use of the formulae for the Wightman 
function and the VEV of the field square, one obtains two equivalent forms, corresponding to 
j = a and j = b (no summation over i): 



A—D 00 



<0|^|0> = (o |^|o ) + m fc ) (j) -^-^^AA 



1=0 



r 

In this formula, for a given function g(z) we use the notations 



F<-%(z)] 
FU[g(z)] 
F®\g{z)] 



H 

1 ( dg(z) 
dz 



3-* 



dg(z) 



dz 



z dz 



+ [ 1 + 

1 

2 



, C d 2 

+ —7-9 ( z ) 

z dz 



UJ 



9\z), 



--5 2 (^) + ^ fi + ^)/r(-.). 



dz 



1 + -0 



5 2 (z) A 2 - m 2 



D - 1 



A? 



(46) 

(47) 
(48) 
(49) 



with (7(2) = zjj\z, Xij), where i = 2, . . . ,D and the indices 0,1 correspond to the coordinates 
t, £, respectively. In formula (|46|1. 



1-D 00 



<Oo|2f |Oo) = tf| 



7T 2 S D 



(50) 



is the corresponding VEV for the vacuum without boundaries, and the term (2^}w) is induced 
by the presence of a single spherical brane located at £ = j. For the brane at £ = a and in the 
region £ > a one has (no summation over i) 



fe\(a) 



1-D 00 



irSn ^— ' 



^AA 2 







ifi (Ajo) 



(51) 
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For the geometry of a single brane at £ = b, the corresponding expression in the region £ < b 
is obtained from (|51(1 by the replacements Q31JI . The expressions for the functions f^[g(z)] in 
(|5Uj) are obtained from the corresponding expressions for i^Wfgfz)] by the replacement to —* ia;. 
It can be easily seen that for a conformally coupled massless scalar the boundary induced part 
in the energy-momentum tensor is traceless. The boundary-free part Q50|) and the single brane 
part (T^)W in the region £ < j are investigated in Ref. [S]. 

Now we turn to the investigation of the brane-induced VEVs in limiting cases. First of all 
let us consider single brane part ()51|) . At large distances from the brane, £ S> rjj, the main 
contribution comes from the I = term and one has 

(T^^X^F^^, (52) 

where ((p 2 )^ is given by ()36|) and 

= 1 - 4C. F« = tf> = 1 - 4 C - i^M. (53) 

In this limit the radial vacuum stress is suppressed by the factor Ao£ with respect to the corre- 
sponding energy density and azimuthal stresses. In the limit a —* when £ is fixed the main 
contribution into the w-integral comes from the lower limit and to the leading order we obtain 

(T?) (a) « - i f 2 /°, : EAA?^[^(A,a=o. (54) 
27r5£) In (r# /a) ^ 

For rij — ► 0, in the way similar to that for the field square it can be seen that for a non- 
minimally coupled scalar field (T^)W is suppressed by the factor exp[— 2yJ C,n{n + !)(£ — a) /Vff]- 
For a minimally coupled scalar the main contribution comes from the I = term and the brane 
induced VEV (EU behaves like r]f D . 

Now let us present the VEV (|46|) in the form 

<0|7f|0> = (Oo|Tf|Oo) + Pi)® + ^) (afe) ' ( 55 ) 

j=a,b 

where the interference part is given by the formula (no summation over i) 



1=0 



iS 6) (A,6)Z w (A,o,A,6) ^ a) (A ; a) 



(56) 



The surface divergences are contained in the single brane parts and the term (|56|) is finite for all 
values a ^ £ ^ b. An equivalent formula for (T^)^ ab ^ is obtained from Eq. Q56|) by replacements 
(j^lj) . The behavior of the interference part (|56|) in the limits a — > and & — > oo is similar to that 
for the field square. In the near-horizon limit, a, b <C r#, for both single brane and interference 
parts replacing the summation by the integration in accordance with formula (|40j) . the result 
for the parallel plates uniformly accelerated through the Fulling-Rindler vacuum is obtained. 



4 Vacuum interaction forces between the branes 

In this section we will consider the vacuum forces acting on the branes. The force acting per 
unit surface of the brane at £ = j is determined by the radial component of the vacuum energy- 
momentum tensor evaluated at this point. By using the decomposition of the VEV for the 
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energy-momentum tensor given by (|46jl . the corresponding effective pressures, p^' = —(T^)^ = j, 
can be presented as the sum 



ij) =P?+P%), j = a,b, (57) 



where the first term on the right is the pressure for a single brane at £ = j when the second 
brane is absent. This term is divergent due to the surface divergences in the subtracted vacuum 
expectation values and needs additional renormalization. This can be done, for example, by 
applying the generalized zeta function technique to the corresponding mode sum. This procedure 
is similar to that used in Ref. [0] for the evaluation of the surface energy for a single brane. 
The second term on the right of Eq. (|57|). P«2t)' * s ^ e P ressure induced by the presence of 
the second brane, and can be termed as an interaction force. This term determines the force 
by which the scalar vacuum acts on the brane due to the modification of the spectrum for the 
zero-point fluctuations by the presence of the second brane. It is finite for all nonzero distances 
between the branes and is not affected by the renormalization procedure. 

For the brane at £ = j the interaction term is due to the third summand on the right of Eq. 
(|4fij). Substituting into this term i = k = l, £ = j and using the Wronskian relation for the 
modified Bessel functions one finds 

l-D oo ,-OQ 

Pgt) = 2^^s~J2 Dl / ^ [( A ^' 2 + ^) % + 4< ^ " !] VjUka, Mb), (58) 
D i=o Jo 

with (3j = Bj/jAj. The interaction force acts on the surface £ = a + for the brane at £ = a and 
on the surface £ = b — for the brane at £ = b. In dependence of the values for the coefficients in 
the boundary conditions, the effective pressures can be either positive or negative, leading 
to repulsive or attractive forces, respectively. For Dirichlet or Neumann boundary conditions on 
both branes the interaction forces are always attractive. For Dirichlet boundary condition on 
one brane and Neumann boundary condition on the other one has Pn^ > and the interaction 
forces are repulsive for all distances between the branes. Note that the interaction forces can 
also be written in another equivalent form 



n U) r l7 D /"» 



J u 1=0 







1+ (4C-1)& 



d , 



£ a) (A,a)^ 6) (W 



lL b) (\ib)K^(X ia ) 



(59) 



Now we turn to the investigation of the interaction forces in the asymptotic regions of the 
parameters. In the limit a — > b the dominant contribution into the expression on the right of 
(|58|) comes from large values I and uj. Replacing the summation over I by the integration in 
accordance with Xy/^o^/(0 ~~ * ^ Jo* dll D ~ 2 f(l)/T(D — 1), and using the uniform asymptotic 
expansions for the modified Bessel functions, to the leading order we find 

n w „ a rmc R (g + i) 

P(int) ~ a *b (47r) (D+l)/2 ( b _ ' ^ DU ^ 

where Cr( 2; ) is the Riemann zeta function, a a b = —1 for 5s a 5s f) = 1 and a a b = 1 — 2~ D for 
^B a ^B b = — 1- Hence, for small distances between the branes the interaction forces are repulsive 
for the Dirichlet boundary condition on one brane and non-Dirichlet boundary condition on the 
other and are attractive for all other cases. Note that in the limit a — > b the interaction part of 
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the total vacuum force acting on the brane diverges, whereas the renormalized single brane parts 
remain finite. From here it follows that at small distances between the branes the interaction 
part dominates. 

When the left brane tends to the horizon, a — > 0, the main contribution into the vacuum 
interaction forces comes from the lower limit of the w-integral and one has 

(a) w s B y H - D f, K^m 

2nS D aHn 3 (r H /a)^ Q Vl jW (A;6 ) ' ^ 

^) ~ 4^Hn 2 (rE/a)^ ' if 2 (A,&) ' 1 j 

In this limit the interaction forces have different signs for the Dirichlet and non-Dirichlet bound- 
ary conditions on the brane £ = a. The combination 0B a P^ t \ is positive for large values j3 b and 
is negative for small values of this parameter. In the limit b — > oo for fixed a, the dominant 
contribution comes from the lowest mode I = and assuming that A b / ±Aoi?b, we have the 
estimates 

(a) Ale-*™ A b -X B b [°° (A 2 a 2 +^ 2 )/3 2 + 4(^-1 

^) ~ 2S D a 2 rx~ 1 A b + ^oBb J *" K^(\ a) ' ^ 



(6) A e- 2A » b Afc - X B b [°° lj a \\ a) 



with the exponentially small interaction forces for both branes. In particular, the combination 
(A* — ^o-^b)^(int) * s positive/negative for large/small values j3 a . For small values of the curvature 
radius r#, in the way similar to that used for the VEV of the field square, it can be seen that 
for a non-minimally coupled scalar field the main contribution comes from the I = term and 
to the leading order we have 

(4\ $R <5r, \Cn(n + l)l 3 ' 4 va& , rz — ; r, % , , s 

pgt) « ~ f^gf 5+175 J- ex P[-VCn(n + l)(6 - a)/r„]. (65) 

For a minimally coupled scalar field the dominant contribution comes from the I = term with 
Ao = m and the interaction forces per unit surface behave like r^j~ D . In the near-horizon limit, 
a, b -C rji, in (|58jl we replace the summation over I by the integration in accordance with (|4()|) 
and to the leading order the result for the geometry of two parallel plates uniformly accelerated 
through the Fulling-Rindler vacuum is obtained. 



5 Conclusion 

In this paper, we investigate the polarization of the scalar vacuum induced by two spherical 
branes in the (D + l)-dimensional bulk Ri x S®^ 1 , assuming that on the branes the field 
obeys the Robin boundary conditions. In the corresponding braneworld scenario based on the 
orbifolded version of the model the coefficients in the boundary conditions are expressed in terms 
of the brane mass parameters by formula © • The most important characteristics of the vacuum 
properties are the expectation values of quantities bilinear in the field operator such as the field 
square and the energy-momentum tensor. As the first step in the investigation of these VEVs 
we evaluate the positive frequency Wightman function. The corresponding mode sum contains 
the summation over the eigenfrequencies. In the region between the branes the latter are the 
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zeros of the bilinear combination of the modified Bessel functions and their derivatives. For 
the summation of the series over these zeros we employ a variant of the generalized Abel-Plana 
formula. This allows us to present the Wightman function as the sum of a single brane and 
second brane induced parts, formulae (|2"T|) and (|27j) . 

The corresponding VEVs of the field square and the energy-momentum tensor are obtained 
from the Wightman function in the coincidence limit and are investigated in section |21 These 
VEVs are given by formulae (|32|) and (|46j) . where j = a,b provide two equivalent representations. 
We have considered various limiting cases of the general formulae. In particular, we have shown 
that when the left brane tends to the horizon the interference parts in the VEVs of the field 
square and the energy- momentum tensor vanish as 1/ In 2 {rjj /a) . In the limit when the right 
brane tends to infinity, b — > oo, the interference parts are suppressed by the factor exp(— 2Ao&). 
In the high curvature regime, corresponding to small values r#, the behavior of the VEVs is 
essentially different for minimally and non-minimally coupled scalar fields. For a non-minimally 
coupled field the VEVs are suppressed by the factor exp[— 2-^/ (n(n + 1)(^ — a)/rn] for single 
brane parts and by exp[— 2-sJ (n(n + l)(b — a)/r#] for the interference parts. For a minimally 
coupled field the main contribution comes from the I = term and the VEVs behave as rlf . In 
the limit when the both branes are near the horizon, to the leading order the VEVs are obtained 
for the geometry of two parallel plates uniformly accelerated through the Fulling- Rindler vacuum. 

In section |1] we have investigated the vacuum forces acting on the branes. These forces 
are presented as the sum of self-action and interaction parts. Due to the well-known surface 
divergences, the self-action part needs additional subtractions. The interaction forces are finite 
for all nonzero interbrane distances and are given by formula (|58|) or equivalently by (|59|) . In 
general, these forces are different for the left and right branes and, in dependence of the values 
of the parameters, they can be either attractive or repulsive. In particular, at small interbrane 
distances they are repulsive for the Dirichlet boundary condition on one brane and non-Dirichlet 
boundary condition on the other, and are attractive for all other cases. When the left brane tends 
to the horizon the interaction forces acting on the left and right branes behave as l/[a 2 ln(rjy/a)] 
and 1/ In 2 (ru/a), respectively. In the limit when the ^-coordinate of the right brane tends to 
infinity, the interaction forces for both branes are suppressed by the factor exp(— 2\ b). 
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